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In this work we show analytically and numerically that wave transport through random waveguides 
can be modeled as a diffusion with an inhomogeneous diffusion coefficient (IDC). In localized regime, 

IDC retains the memory of the source position. In an absorbing random medium, IDC becomes 
independent of the source. 

PACS numbers: 42.25.Dd,42.25.Bs,72.15.Rn 


I. INTRODUCTION 

Diffusive description of wave transport random 
medium has a long historyi. This macroscopic approach 
describes the ensemble-averaged intensity of the wave on 
scales longer than transport mean free path i. As such, 
it has a tremendous practical advantage compared to the 
direct solution of wave equation for each statistical re¬ 
alization of disorder and subsequent averaging over the 
ensemble of solutions. 

Diffusion coefficient can become renormalized^ due to 
wave localization phenomenon^. In three-dimensions for 
sufficiently strong disorder, the diffusion vanishes for 
an infinitely large system^. In practice, however, one 
deals with transport through finite systems. In pro¬ 
cess of adapting self-consistent theory (SCT) of local¬ 
ization^!^ and super-symmetric (SUSY) theories^, it be¬ 
came cleari2iS that the diffusive-like description can also 
be applied to the finite systems that exhibit localized 
transport, in particular, to the low dimensional systems. 
In such description the diffusion coefficient becomes in¬ 
homogeneous in space and dependent on system size^r^i 
and geometryi^. In quasi-lD or ID lossless medium both 
SCT and SUSY lead to the following equation for the 
ensemble-averaged intensity (/ {z,z')) in the presence of 
a point source Jq at z' 

— {d/dz)D{z)d{I {z, z'f) jdz = (z — z') (1) 

Diffusion of this kind results leads to a highly unusual 
macroscopic transporlji^. We stress that the medium it¬ 
self, i.e. the density of scatterers, is uniform and that the 
inhomogeneous diffusion is brought about by non-local 
wave interference effects. 

To date, the studies of inhomogeneous wave diffusion 
has concentrated on geometry where a wave is incident 
onto the random medium from an outside (free space) 
regioit^ ~^^d^~ — . Inhomegeneous wave diffusion descrip¬ 
tion was successful in describing the light intensity under 
various measurement conditions, see Ref.— for a review. 
Ensemble-averaged intensity, however, is only the first 
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step in characterizing wave transport in random media. 
Indeed, the second order statistical quantities, such as 
fluctuations or correlations, become important at the on¬ 
set of Anderson localizatio n^^d^ , they require the knowl¬ 
edge of the green function of the diffusion equation, e.g. 
(I (z, z')) with arbitrary z—^—. 

In this work we test the applicability of Eq. [T] with 
an arbitrary position of the source. We show analyti¬ 
cally that for z' ^ 0, the diffusion equation is applicable, 
however inhomogeneous diffusion coefficient (IDC) D{z) 
acquires the dependence on the position of the source z'. 
We derive an closed form analytic expression for D [z, z') 
and verify it with ab-initio numerical simulations. We 
show that D (z, z') is reduced to the known resuHj^i for 
z' —^ 0, i.e. for wave incident from the outside region. 
We demonstrate that when absorption, unavoidable in 
an experiment, is present in the system, D[z,z') looses 
its dependence on source position and that it can be de¬ 
termined using self-consistent theorySiiS. 

This paper is organized as follows. In Sec. |n] we ob¬ 
tain an analytical expression for IDC that describes wave 
transport in a single mode waveguide with an external 
source. We verify the applicability of the result with nu¬ 
merical simulations. In Sec. uni we demonstrate analyt¬ 
ically and confirm numerically that IDC depends on the 
position of the source inside random medium. The effect 
of absorption is accounted for in Sec. El where we show 
that IDC looses is position dependence for a sufficiently 
strong absorption in the system. 

II. INHOMOGENEOUS DIFFUSION IN 
PASSIVE SYSTEM WITH AN EXTERNAL 
SOURCE 

We consider one-dimensional random medium occupy¬ 
ing Q < z < L region. Propagation of a scalar wave E{z) 
is described by 

l^ + k^[l + e{z)]E{z) = Q, (2) 

here k = 27r/A is wave number and e{z) is a random 
process. For a wave with unit amplitude incident from 
the left, the boundary conditions can be expressed in 
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terms of reflection r and transmission t coefficients as 

E{z) = + r z < 0 

E{z) = z>L (3) 

We are interested in obtaining a closed form expression 
for the intensity (/(z)) = (|£l(z)|^) averaged over the 
ensemble of random processes e(z). Indeed, inhomoge¬ 
neous diffusion coefficient D(z) can be found from the 
Pick’s law 

(J(z)) = -D{z)d{Iiz))/dz, (4) 

where (J(z)) is the flux. In passive random medium the 
flux is conserved during propagation and, thus, can be 
found from the boundary conditions. Indeed, the fraction 
of the flux propagating in positive (-I-) or negative (—) 
direction can be expressed asi 

(J(±)(z)) = v{I{z)) T {D{z)/ 2) d{I{z))/dz (5) 


where v is the wave speed. The right boundary (z)) 
vanishes so {J{z)) = {J^^\L)) = 2v{I{L)). Substituting 
this expression into Eq. o we obtain 

D{z) = -2v{I{L))/{d{I{z))/dz). (6) 

Therefore, finding IDC requires the knowledge of (only) 
(/(z)) for the problem defined by Eqs. (I2I3I) . Such so¬ 
lution has been obtained in Refsj^^— . The common 
theme is these studies is to relate the statistical prop¬ 
erty of the wave field inside the medium to those at the 
boundary where it can be obtained using the limit¬ 
ing theorems. Such approach is in spirit of the well 
known self-embedding method^^. We will assume that 
e(z) is a delta-correlated Gaussian process with (e(z)) 
and (e(z)e(z')) = a5{z — z'). Under these conditions the 
solution for the ensemble-averaged intensity is obtained 
in the fornt2^“— 


(7(z)) = I- 



iC-iz-L/2)/0^ - 

L/^ 




dc. 


( 7 ) 


where we introduced localization length as = afc^/2. 
Substitution of Eq. JT]) into Eq. ([51) gives us the analytical 
expression for IDC. 

A compact expressions for both (/(z)) and D{z) can 
be obtained when L ^ In this limit, the expression in 
parentheses in the integrand of Eq. © can be approxi¬ 
mated with a step function h{^) and the integral can be 
computed in terms of the error function erf(x): 

(/(z)) ~ I-erf [5c] (8) 

with scaling parameter Zc = (z — L/2)/i/Lj as the argu¬ 
ment. In Fig. [T^, (/(z)) is plotted for L/^ = 5, 10, 20, 
50 and 100 with and without scaling z-coordinate. We 
confirm that Eq. (I5|) approximates the exact expression 
Eq. ([71) well. We note, that such distribution has been 
observed in numerical simulations of energy deposition 
in wavefront shaping in random medium^l. 

Asymptotic expression IDC in L ^ is obtained by 
substituting Eq. ([5|) into Eq. (jH]): 

D(z) ~ Do exp [5^ - L/4^] . (9) 

Here Do = vi is unrenormalized value of diffusion coeffi¬ 
cient in terms of transport mean free path 7 = 2^. Fig.[TjD 
confirms the universality of IDC inside ID passive ran¬ 
dom medium in terms of the scaling parameter Zc- Eq. 
agrees with that derived in super-symmetric theory of 
RefJ^ for quasi-lD geometry (a multimode waveguide). 


III. INHOMOGENEOUS DIFFUSION IN 
PASSIVE SYSTEM WITH AN INTERNAL 
SOURCE 

For the source located inside random medium, Eq. © 
is modified to include a point source at z'\ 

^ ^ ^ + e(z)] E{z, z') = 5{z - z'), (10) 

whereas the boundary conditions in Eq. m are replaced 
with the outgoing wave conditions at both ends of the 
waveguide 

E(z) = ti z < 0 

E{z) = t 2 z>L (11) 

Under these conditions, the flux inside the medium is 
a piece-wise constant function with a jump at the posi¬ 
tion of the source z'. The values of (J(z, z')) for z < z' 
and z > z' can be determined by applying Eq. © at 
z = 0 and z = L respectively. We find (J(z < z')) = 
—2u(/(0,z')) and (J(z > z')) = 2v{I{L,z')), where 
(/(z, z')) = (|U(z,z')|^) and E{z,z') is the solution of 
Eq. (11011 with boundary conditions in Eq. (HU). There¬ 
fore, IDC can be written based on Pick’s law Eq. as 

= -2.(d(7(.))/d.)-' I ^ 

( 12 ) 
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As in Sec. im the above expression for IDC requires the 
knowledge of the ensemble-averaged intensity {I{z,z')). 
The latter has been obtained in Refs.— >2^ in the form 



FIG. 1: (a) Spatial distribution of ensemble-averaged inten¬ 
sity {I{z)) from Eq. (O (blue lines) without scaling (inset) 
for systems with L/^ — 5,10,20,50,100 (left to right). The 
main panel presents the same data in terms of scaled co¬ 
ordinate Zc = {z — L12) j the approximate expression 
given by Eq. ((S]) is shown with thick dashed line, (b) Inhomo¬ 
geneous diffusion coefficient found from Eq. m for systems 
with L/^ = 5,10, 20, 50 without (inset) and with (main panel) 
scaling. Thick dashed line is found from asymptotic L ^ 
expression in Eq. ©. 


J 


(/(z, z')) = TT exp 


3L 


d/i exp 




sinh(7ru) ,, , 1. , (z + z' — L). 

-^ + 7 cos --- ’-) 

/iCOsh(7r/x)^ 4 ^ 


2 1. , (L—\z — z'\). , . (L—|z —z'l), 

+ {fj, —-)cos(2/r--- )+fisin{2fi ---^ 


(13) 


Substituting this expression into Eq. (Ha gives us the 
final result. 

We make the following observations. First of all, in 
the limit of z' = 0, we recover the result for an external 
source found in the previous section. Indeed, Eq. CSD can 
be showit^ to reduce to Eq. 0. Second, unlike Eq. o, 
IDC D(z, z') depends on the source position z'. In Fig. [2] 
we evaluated Eq. (na) for L/^ = 7.6 and four values of 


z': 0 (outside source), i/4, i/2 and 3i/4. Indeed, IDC 
shows strong dependence on the position of the source. 
We note that D{z, 0) is always greater than D{z, z' > 0) 
with the minimum value at the middle of the sample for 
z' = i/2. 

We verified the above results with numerical simu¬ 
lations of wave with k = 1.45 propagating normally 
through a stack of alternating dielectric slabs with di- 
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FIG. 2: (Color online) Inhomogeneous diffusion coefficient 
inside ID passive random media of L = 7.6^ computed for 
four source positions z' = 0 (upper curves), I//4, and 
fiLjA:. The latter three curves have cusps at the position of the 
source. Dashed lines were found by substituting the analytical 
result Eq. in into Eq. (HU). The solid lines were obtained 
numerically. 


electric constants ei = 1 and e 2 = 1.2. The widths 
of the stacks of the first kind is distributed uniformly 
in the interval di G (0.9,1.1), while the width of the 
others is being kept constant d 2 = 1- The propaga¬ 
tion of the waves in the system consists of free prop¬ 
agation in the slabs and scattering at the interfaces, 
where the boundary conditions should be satisfied. It 
can be described using the transfer-matrix formalism, c.f. 
Ref i^9i30 . We computed J(z) = —kcIm[B(z)dB(z)/dz] 
and I(z) = (k^/2) \E{z)\^ -f (1/2) \dE{z)/dz\^ in a sys¬ 
tem with N = 8 X 10^ layers numerically and averaged 
over 10® disorder realizations. IDC was found from Pick’s 
law Eq. ®. The results reported in Fig. [5] (sold lines) 
agree with the analytical expression (dashed lines). 


IV. EFFECT OF ABSORPTION 

Absorption is inevitable in optical experiments. The 
effect of absorption is to suppress resonant tunneling of 
the waveii, thereby to increase In this section 

we perform numerical analysis of the effect of absorption 
on IDC for an internal source. 

Modeling absorbing random medium is accomplished 
by adding a constant imaginary part to e{z) in Eq. OHD 
as e{z) e{z) + iy. Addition of the loss results in 
an extra term ^~^{I{z,z')) in the inhomogeneous dif¬ 
fusion equation. Absorption length can be obtained 
for a given value of 7 from the continuity condition 
d{J{z,z'))/dz = {I{z,z'))/Ta where Ta 

FigureOshows IDC obtained numerically for the model 
in Sec. uni We choose number of layers in a stack N = 



zjL 

FIG. 3: (Color online) Inhomogeneous diffusion coefficient 
in absorbing random medium with = 15.2. The lower 
set of curves corresponds to absorption = 3.2 and four 
positions of the source. Thick red line is the prediction of 
SCT from Eq. (I14I15II . For stronger absorption Lj^a = 10 
(upper set), SCT predicts IDC well. 


1.6 X 10"^ and two values of 7 : 10“® and 10“^. These pa¬ 
rameters give Lj^ao- 3.2 and 10 respectively. Similar to 
the passive (non-absorbing) case in Sec. IIIII IDC clearly 
shows dependence on the position of the source inside 
the medium {z' = 0, L/4, T/2 and 3L/4 are shown), 
it has a cusp feature at z'. However, closer inspection 
of Fig. [3] shows that dependence on source position is 
strongly suppressed at large T/^a; in this limit be¬ 
comes comparable to the localization length 

Performing computationally expensive numerical sim¬ 
ulations is not practical, in particular, in higher dimen¬ 
sional systems. Self-consistent theoryiiSiiiS has been suc¬ 
cessful in providing a good prediction for IDC for systems 
with L/^ not too largeii, and it was shown^ to be ac¬ 
curate in the absorbing systems. In all previous works, 
external source has been considered. Here, we computed 
the prediction of the self-consistent theory for systems 
with different amount of absorption. Self-consistency 
condition relates IDC to the return probability {I{z,z)) 
as 

Do/D(^^^Hz) = 1 + v{I{z,z)). (15) 


These equations form a closed set, sufficient for finding 
jjiSCT)^^y The solution of Eqs. (I14I15I) is shown with 
a thick solid line in Fig. |31 We find that in presence of 
sufficiently strong absorption SCT makes adequate pre¬ 
diction for IDC even for internal sources. 
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V. CONCLUSION 

In this work we investigated the applicability of the 
inhomogeneous diffusion approach to describing wave 
transport in random medium. We have shown analyt¬ 
ically and numerically that in the regime of localized 
transport the inhomogeneous diffusion coefficient 

can be defined through Pick’s law. The benefit of such 
approach, is that it allows one to obtain the ensemble- 
averaged value of intensity without the need to perform 
statistical averaging. 

Our analysis shows that IDC exhibits a significant de¬ 
pendence on the source position z'. Such dependence 
has not been discussed before. This is because previous 
studies have concentrated on the common experimental 
arrangement - the incident wave is impingent onto the 
sample from the outside. Present study of the IDC with 
an internal source is of practical interest for a number of 
reasons. First of all, solution {I{z,z')) of the diffusion 
equation with IDC with an internal source z' is the green 


function which can be used to define second order statis¬ 
tics (e.g. fluctuation, correlations) of wave transport. 
Secondly, Pick’s law with IDC in the form of D{z,z') 
points to a highly unconventional type of diffusion in the 
localized systems. The spatial dependence of IDC has 
been showni^ to exhibit an unusual macroscopic trans¬ 
port behavior. The additional dependence on the source 
position found in this work may necessitate completely 
new non-local approach to transport. Accurate descrip¬ 
tion of wave transport through random media would in¬ 
form the studies of limits of wavefront shaping^ with 
applications, in particular, in the field of biological imag- 
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